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1^^ ■ Abstract It is well-known that the celebrated Camassa-Holm equation has the peaked 

CA ■ solitary waves, which have been not reported for other mainstream models of shallow 

water waves. In this letter, the closed-form solutions of peaked solitary waves of 
lAl \ the KdV equation, the BBM equation and the Boussinesq equation are given for the 

first time. All of them have either a peakon or an anti-peakon. Each of them exactly 
2 ' satisfies the corresponding Rankine-Hogoniot jump condition and should be understood 

as weak solution. Therefore, the peaked solitary waves might be common for most of 

shallow water wave models, no matter whether or not they are integrable and/or admit 
•r^ ' breaking-wave solutions. 
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Since the solitary surface wave was discovered by John Scott Russell [13] in 1834, 

many models of solitary waves in shallow water have been developed, such as the 
K> I Boussinesq equation ^, the Korteweg & de Vries (KdV) equation p^, the Benjamin- 

H ■ Bona-Mahony (BBM) equation [T], the Camassa-Holm (CH) equation |1], and so on. 

■ ■ ■ Unlike other models, the celebrated CH equation 

Ut + 2kUx — Uxxt + 'iuUx = 2UxUxx + UUxxx (1) 

possesses both of the global solutions in time and the wave-breaking solutions, where 
u denotes wave elevation, x, t are the spatial and temporal variables, the subscript 
denotes the differentiation, and < k < 1/2 is a constant related to the critical 
shallow water wave speed, respectively. Especially, when k = 0, the CH equation ([1]) 
has a peaked solitary wave u{x,t) = c exp(— |x — c t|), where c denotes the phase 
speed. It should be emphasized that, at the wave crest x = at, the above closed-form 
solution has no continuous derivatives respect to x, thus it does not satisfy Eq. ([1]) 
(when K = 0) at the crest x = at. So, it is not a strong solution. However, Constantin 
and Molinet |6] pointed out that u{x,t) = c exp(— |s — c t|) could be understood as 
a weak solution of the CH equation ([1]) when k = 0. 

Note that such kind of discontinuity (or singularity) exists widely in natural phe- 
nomena, such as dam breaking [IT] in hydrodynamics, shock waves in aerodynamics. 



black holes described by the general relativity, and so on. In the frame of water waves, 
Stokes [H] found that the limiting gravity wave has a sharp corner at the crest. In 
case of the dam breaking, there exist sharp corners of wave elevation at the beginning 
t = 0, and such kind of discontinuity of the derivative of wave elevation does not 
disappear due to the neglect of the viscosity, as shown by Wu and Cheung [15]. In 
fact, problems related to such kind of discontinuity belong to the so-called Riemann 
problem |2lll2l[T5] , which is a classic field of fluid mechanics. 

However, such kind of peaked solitary waves have never been found for other 
mainstream models of shallow water waves. Solitary wave solutions in shallow water 
with a fixed speed c and permanent form were found by Korteweg & de Vries ^U] 
using the KdV equation 

Ct + Cxxx + 6C Cx = 0, (2) 

subject to the boundary conditions 

C -> 0, Cx -> 0, Cxx -^ 0, as |x| -> +00, (3) 

where C, denotes the wave elevation. The KdV equation ([2]) admits solitary waves but 
breaking ones. The traditional solitary wave of KdV equation (|2]) reads 



c .2 



C{x,t) = - sech 



— {x-ct-xo) 



(4) 



where sech denotes a hyperbolic secant function and Xq may be any a constant. This 
solitary wave has a smooth crest with always positive elevation ({x, t) > 0, and besides 
its phase speed c is dependent upon the wave height, i.e. c = 2(rnax, so that higher 
sohtary waves propagate faster. To the best of the author's knowledge, no peaked 
solitary waves of KdV equation (|2]) were reported. 

The KdV equation ([2]) is an approximation of the fully nonlinear wave equations 
that admit the limiting gravity wave with a peaked crest, as pointed out by Stokes [Ti] . 
Besides, as mentioned before, this kind of discontinuity of wave elevation widely exists 
in hydrodynamic problems such as dam break pT]. In addition, current investigations 
[5l[7] reveal the close relationships between the CH equation ([1]) and the KdV equation 
([2]). Therefore, one has many reasons to assume that the progressive solitary waves 
of the KdV equation ([2]) might admit peaked solitary waves, too. 



2 Peaked solitary waves of KdV equation 

Write ^ = x — ct — Xq and ri{^) = ({x,t). The original KdV equation ([2]) becomes 

— CT]' + T]"' + Grjr]' = 0, (5) 

subject to the boundary condition 

V^O, r/' ^ 0, 7]" -^ 0, as 1^1 -^ +00, (6) 



where the prime denotes the differentiation with respect to ^. Besides, there exists 
the symmetry condition 



v{-0 = v{0, ee(-oo,+oo) 



(7) 



so that we need only consider the solution in the interval ^ > 0. Integrating (|5]) with 
(E]) gives 



- cr] + r]" + 3r]^ = 0, ^ > 0. 



Multiplying it by 2r]' and then integrating it with ([6]), we have 
which has real solutions only when 



< 



Under the restriction flTO|) . we have from (|9]) that 



say, 



r]' = ±r] a/c — 2r] 
drj 



7] a/c-2?7 



±d^, e > 0. 



(8) 
(9) 

(10) 

(11) 
(12) 



Let us first consider the solitary waves with f] > 0. Integrating (fT2l) in case of 
< ?7 < c/2, we have 



tanh 



2r] 



±^^ + «, 



(13) 



where a is a constant. Since < 21] /c < 1, the left-hand side of the above equation 
is non-negative so that it holds for all ^ > if and only if 



tanh 



1_^ 
c 



^e + «, e>0,a>0. 



which gives 



^(0=2 s^^^^ 



2 ^ 



Using the symmetry condition (^^, it holds 



^(0 = 2 ^^^^' 



■iei + « 



, e>0,a>0. 



-oo < ^ < +c>o, a > 



in the whole interval. Thus, we have the solitary waves of the first kind 



C(x,t) = -sech^ 



— \x — c t — Xo\ + a 



(14) 



(15) 



(16) 



where a > 0, with the wave height 



[a), a > 0. 



:i7) 



Note that a > is a constant parameter. When a = 0, it is exactly the same as the 
traditional solitary waves of the KdV equation with a smooth crest. 

In case of r^ < 0, it holds 1 — 21] /c > 1 so that 

tanh~ 



1_^ 
c 



becomes a complex number which has no physical meanings. This is the reason why 
the solitary solutions of the KdV equation in case of ?7 < was traditionally neglected. 
So, we must be very careful in this case. 



In case of r^ < 0, write 



v/c-2r/ = y^./l-^ 



'c z. 



- = i-H!i>i. 



T] = -(1 — z^), drj = —c z dz. 



where z > 1 and 
Then, we have 

Since 2; > 1, it holds 

dr] 
risjc — If] ^\z — \ -2 + 1 

Substituting it into ( TT^ and integrating, we have 

' z-\ 



dz dz 






In 



^ + 1 



±v^e-2/3, e>o. 



where /3 is a constant. Since 2; > 1, it holds 



which gives 



0<^<l. 

2; + 1 



In ( 1 < 0. 



^ + 1 



Thus, ([20]) holds for all ^ > if and only if 

' z-\ 



In 



z+1 



ce + 2/3), e>0, /3>0, 



(19) 



(20) 



(21) 



which gives 



coth 



'V~c 



e + /3 



^> 0, 13 >Q. 



So, we have 



^(0 = 2(1-^') 



— ( - ) csch^ 



'^i^^ 



e>0, /3>0, 



where csch denotes a hyperbohc cosecant function, and /3 > is a constant. Using 
the symmetry condition ([7]), we have the sohtary wave 



Vii) = 7t(1 



- I csch^ 
2. 



^ki + /j 



-cx) < ^ < +00, /3 > 0, 



which is vahd in the whole interval. 

Let d < ({x,t) < denote the restriction of ({x,t). Then, it holds 

viO) = - (I) csch2(/3) > d, 

which gives 

2d 



(5 > csch-M . 

Thus, we have the solitary waves of the second kind: 



C(x,t) = --csch^ 


-^^ a; ■ 
. 2 


-ct- 


-a;o +/3 


where 

r 






13 > csch-M- 


2d 
c 




with the restriction d < C{x,t) < 0. 









(22) 



The closed- form solutions fll6p and (B2]) of the solitary waves of the first and second 
kind exactly satisfy the KdV equation ([2]) in the whole domain —00 < a; < +00 and 
t > 0, but except X = ct + xo (corresponding to the wave crest). This is rather similar 
to the closed-form solution u{x,t) = cexp(— |a: — c t|) of the CH equation ([1]) when 
K = 0. Thus, the closed-form solutions f lT6|) and fl22|) should be understood as weak 
solutions. Besides, the corresponding Rankine-Hogoniot jump condition pJJ must be 
satisfied, so as to ensure that f lTBl) and fl22]) have physical meanings. 



To give the corresponding Rankine-Hogoniot jump condition [TTj of the KdV 
equation ([2]), we rewrite it in the flow 



Ct + [/(C)]. = 0, 

where f{() = (^x + 3C^- Then, C is a weak solution of 



if 



+00 



[(Vt + fiOVx] dxdt + 



C(a;, 0)(^(a;, 0)dx = 



(23) 



(24) 



for all smooth functions if with compact support. Besides, there exists such a theorem 
that, if (^ is a weak solution of ( 123|) such that C is discontinuous across the curve 
X = a{t) but (^ is smooth on either side oi x = a{t)^ then C, must satisfy the condition 

i/±\ J \^ ) J \^ ) /OCX 

^ w = — (- -c+ — ' 

across the curve of discontinuity, where (~{x,t) is the limit of ( approaching {x,t) 
from the left and (~^{x,t) is the limit of ( approaching {x,t) from the right. 

Then, due to the symmetry of the progressive peaked solitary waves about the 
crest X = c t + xo, we have 

/(r ) - /(c+) icu - icu + 3(c-)' - 3(C)^ 



C--C+ C--C 



C--C+ 

(,S> )xxx [^ 



(C-). - (C+) 



- + 3 (C" + C^) 

+ 3 (C" + C+) , (26) 



which provides us the so-called Rankine-Hogoniot jump condition of the KdV equa- 
tion: 

l,^ Jxxx Is J 



C 



(C-). - (C 



+ 



+ 3 (C" + C"^) , as X -> c t + Xq. (27) 



/x 



For the closed-form solution flT6ll at the crest, we have 



C v2, 



r = C = -sech\a), (28) 

3/2 

(C")x = ^sech2(«)tanh(a), ((+), = -((-),, (29) 






iOxxx = — [cosh(2a) -5]sech (a)tanh(a), {C^)xxx = -{C)xxx- (30) 

Substituting them into fl27|) . we have 

- [sech^ia) + sech2(a)cosh(2a) - 2] = 0, (31) 

which is an identity for not only a = (corresponding to the traditional smooth 
solitary wave) but also arbitrary a > (to the peaked solitary wave)! 

Similarly, for the closed- form solution f l22l) at the crest, we have 

(32) 
(33) 

:• = -iOxxx. (34) 



C" = 


- C+ = -^csch2(/3), 


(C-). = 


„3/2 

= __csch2(/3)coth(/3), (C+). = -(r)., 


Is )xxx - 


p5/2 

[cosh(2/3) + 5]csch^(/3)coth(/3), (C+ 




X - c t - x„ 



Figure 1: The solitary waves C{x,t) of the first kind of the KdV equation (j2]) with 
the same phase speed c = 2. Red hne: a = 0; Green hne: a = 1/2; Blue line: a = 1; 
Black line: a = 3/2. 




X - c t - x„ 



Figure 2: The solitary waves of the second kind of the KdV equation ([2]) with the 
same phase speed c = 2. Red line: /3 = 1; Green hne: (3 = 5/4; Blue line: /3 = 3/2; 
Black hne: (3 = 7/4. 



Substituting them into ( !27l) . we have 



^ [csch\(3) cosh(2/3) - csch^{f3) - 2] = 0, 



which is an identity for arbitrary constant /3 > 0. 



(35) 



Therefore, both of (11611 and (1221) satisfy the corresponding Rankine-Hogoniot jump 
condition and should be understood as weak solutions of the KdV equation (|2]). It 
should be emphasized that the traditional solitary wave (jl]) with a smooth crest is just 
only a special case of the solitary waves flT6|) of the first kind when a = 0. Note that 
the solitary waves fITBl) and fl2^ of the first and second kind of the KdV equation ([2]) 
have a peakon and an anti-peakon, respectively. Since a > and P > are arbitrary 
constant, the phase speed of these peaked solitary waves has nothing to do with the 
wave amplitude, for example as shown in FigsiT]and[2l This is quite different from the 
traditional smooth solitary wave. To the best of the author's knowledge, such kind of 
peaked solitary waves with peakon or anti-peakon have never been reported for the 
KdV equation. All of these reveal the novelty of the peaked solitary waves (IT6l) and 



m. 



3 Peaked solitary waves of modified KdV equation 

Similarly, the modified KdV equation [HI 



Ct + Cxxx ± GCCx = (36) 

has the two kinds of solitary waves 

2c 

C(a;,t) = ±^ ; ^ ; — , (37) 

^ ' p\/c\x—ct—xo\+a^Qp—\/c\x—ct—xo\—a ^ ' 

where c is the phase speed and a > is a constant. 

Note that ([25]) holds for ([23]) in general. Now, we have /(() = Cxx±2C^. And the 
corresponding Rankine-Hogoniot jump condition reads 

/(C-) - /(n 



C-- 


-{C)xx±2[{Cf-{Cf] 


Is )xx 


C--C+ 



s ~ s 
which provides us the Rankine-Hogoniot jump condition of the modified KdV equation 

(,S> )xxx \S )i 



(C-). - (C+). 



± 2 [((-)' + C"C+ + (C+)'] , as X ^ c t + xo. (39) 



Using the closed-form solution (137)) . it is found that the above condition is satisfied 
for arbitrary constant a > 0. Therefore, the closed-form solution ( 137)) exactly satisfies 
the corresponding Rankine-Hogoniot jump condition of the modified KdV equation 
fl36l) and should be understood as a weak solution. 



4 Peaked solitary waves of the BBM equation 

Similarly, the BBM equation [T] 

Ct + Cx + CCx- Cxxt = (40) 

has the solitary waves of the first kind 



C{x,t) = 3(c-l)sech2 



v^ 



c 



-1 



iei + « 



a>0, (41) 



and the solitary waves of the second kind 



C{x,t) = — 3(c— l)csch^ 



VT 



^-1 



iei + /3 



/3>0, 



(42) 



where ^ = x — ct ~ Xq, cis the phase speed, and a > 0, (3 > are constants. 

Now, we have /(C) = C + C^/2 — Cxt- And the corresponding Rankine-Hogoniot 
jump condition reads 

/(C-)-/(C+) 

^-c+ 
r ) - (c+) + [(c-)^ - icr] /2 - [(r ).* - (o..] 



C--C+ 

2^^ +^ ^ (C-).-(C+). ^ 



which provides us the Rankine-Hogoniot jump condition of the MMB equation 



_ 1 , _ (/■- I /■+\ _ iC )xxt (,C jxxt 

2^^ +^ > (C-).-(C+). ^ 



as X — 7- c t + xq. 



(43) 



(44) 



For the sohtary waves (HTj) of the first kind at the crest, we have 

C~ = 3(c — l)sech^(a), 
((-)^ = 3c(l-c"^)^/^sech2(a) tanh(a), 

iOxxt = -^c\l - c-^f^^sech\a) [sinh(3a) - llsinh(a)] 



(45) 
(46) 

(47) 



and 



C — C ) [C )x — — (C )x, (C )xxt — — (C )xxt- 

Substituting all of these into f j44l) . it is found that ( l4T|) satisfies the Rankine-Hogoniot 
jump condition (jHj) for arbitrary constant a > 0. Similarly, ( H2l) satisfies the Rankine- 
Hogoniot jump condition ( l44l) for arbitrary constant /3 > 0. Therefore, both of ( 14T1) 
and fH2|) should be understood as weak solutions of the BBM equation pUj) . 



5 Peaked solitary waves of Boussinesq equation 



In addition, the Boussinesq equation [3] 



^ -ghTT^-gh 



dt'^ " dx^ 
has the solitary waves of the first kind 



3C^ /^5^ 
dx^ \ 2h 3 dx"^ 



(4^ 



C = /i ( ^ - 1 ] sech^ 



2h\] gh 



-1 l^l + a 



(49) 
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and the solitary waves of the second kind 
C = -h (^-1 I csch^ 



^Vs^-^i^i"^ 



(50) 



where ^ = x — ct — xq, h denotes the water depth, g the acceleration of gravity, c the 
phase speed of wave, a > 0, (3 > are constant, respectively. 

For progressive solitary wave, we have 

d__ _ld_ 
dx c dt 

Then, fl48l) can be rewritten in the form 



d^C gh d^C gh d^ /3C^ h^ d^C 
df^ c dxdt c dxdt \ 2h 3 dx'^ 



0. 



(51) 



Integrating it with respect to t and using the boundary condition (^(±oo) = 0, we 
have 



> gh 

Ct + — 

c 



2h 



(52) 



Thus, we have here 



/«'-?(c-f4<-)- 



According to fl25l) . we have 



C--C 
^^1 (r-n + i[(r)^-(c+)^] + ¥ [(c 



(C^ 



gh 

c 
gh 

c 



^+2X(^" + ^^) 



C--C- 

h^\ {CU - (C+). 

3 



C--C- 



i + |^(r + c+) + fP ^^'^^'^'^^ ^^""^^ 



(53) 



3 J (c-)x. - (C+). 

which provides us the Rankine-Hogoniot jump condition of the Boussinesq equation 



^-l + ^(C- + C) + (y^ (C-).-(C- 



/,2\ (/■-) _ (/■+) 

"' \ VS Jxxx \S /a 



as X — )> c t + Xq. (54) 



For the peaked solitary wave (1491) . we have at the crest that 



(c- 



hi— — 1 I sech^(a), 

v3 I ^ — 1 ) sech^(a) tanh(a) 



(c- 



3^3 f c 



^-1 



4/2,2 ygf^ 



5/2 



sech^(a) [sinh(3a) — llsinh(a)] 



(55) 
(56) 

(57) 
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and 

Substituting all of these into fl5^ . we have 



- I — — 1 I Tsech^fa;) + sech^(a) csch(a) sinh(3a) — 4] = 0, 
4 \gh J 



(5^ 



which is an identity for arbitrary constant a > 0. Similarly, it is found that fl50|) 
satisfies the Rankine-Hogoniot jump condition flMl) for arbitrary constant /3 > 0. 
Thus, the peaked solitary waves fH^ and fISUl) should be understood as weak solutions 
of Boussinesq equation fj48|) . 



6 Conclusion and discussion 

In this article, we give, for the first time, the closed-form solutions of the peaked 
solitary waves of the KdV equation [TU], the modified KdV equation [SlIIS], the BBM 
equation [1], and Boussinesq equation [3], respectively. All of them exactly satisfy 
the corresponding Rankine-Hogoniot jump condition for arbitrary constant a > 
or /3 > 0. Note that the Ist-derivative C,x of the elevation has a jump at the crest. 
Obviously, for all peaked solitary waves found in this article, (^a, > on the left of the 
crest, but Cx < on the right of the crest, respectively. In other words, C,~ > Ct ■ ^O; 
these peaked solutions also satisfy the so-called entropy condition. Thus, they could 
be understood as weak solutions. 

Note that the traditional smooth solitary waves are just special cases of the peaked 
solitary waves when a = 0. So, these weak solutions are more general. Besides, unlike 
the smooth solitary waves whose phase speed strongly depends upon wave amplitude, 
the peaked solitary waves (when a > or /3 > 0) have nothing to do with the wave 
amplitude, for example as shown in Figsil] and [2l This is quite different from the 
traditional solitary waves. In addition, the solitary waves with negative elevation 
such as shown in Fig. [2] have never been reported for these mainstream models of 
shallow waves. All of these show the novelty of these peaked solitary waves. If these 
peaked waves as weak solutions indeed exist and have physical meanings, then nearly 
all mainstream models for shallow water waves, including the KdV equation [10], the 
modified KdV equation [8l[T6], the BBM equation [1], and Boussinesq equation [3] 
and the CH equation [3], have the peaked solitary waves, no matter whether or not 
they are integral and admit breaking-wave solutions. Thus, the peaked solitary waves 
might be a common property of shallow water wave models. As shown currently 
by Liao [9], even the exactly nonlinear water wave equations also admit the peaked 
solitary waves: this might reveal the origin of the peaked solitary waves in shallow 
water predicted by these mainstream models. 

Certainly, further investigations on these peaked solitary waves are needed, espe- 
cially the stability of them, the interactions between multiple peaked solitary waves, 
and so on. Note that these peaked solitary waves have discontinuous 1st derivative 
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at crest so that their higher derivatives at crest tend to infinity and the perturbation 
theory does not work. Thus, strictly speaking, the vahdity of the KdV equation [10] , 
the modified KdV equation p[T6]. the BBM equation [Ij, and Boussinesq equation [3] 
should be verified carefully in the meaning of the weak solution ( l24l) . 



Note that the peaked solitary waves have never been observed experimentally and 
in practice. Obviously, it is more difficult to create and remain such kinds of peaked 
solitary waves than the traditional smooth ones. So, it is an challenging work to 
observe these peaked solitary waves in experiments. 

It should be emphasized that the discontinuity and/or singularity exist widely in 
natural phenomena, such as dam break in hydrodynamics, shock waves in aerodynam- 
ics, black holes described by general relativity and so on. Indeed, the discontinuity 
and/or singularity are rather difficult to handle by traditional methods. But, the 
discontinuity and/or singularity can greatly enrich and deepen our understandings 
about the real world, and therefore should not be evaded easily. 
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